Topological Phase Transition and Electrically Tlinable Diamagnetism in Silicene 
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Silicene is a monolayer of silicon atoms forming a honeycomb lattice. The lattice is actually made of two 
sublattices with a tiny separation. Silicene is a topological insulator, which is characterized by a full insulating 
gap in the bulk and helical gapless edges. It undergoes a phase transition from a topological insulator to a 
band insulator by applying external electric field. Analyzing the spin Chern number based on the effective 
Dirac theory, we find their origin to be a pseudospin meron in the momentum space. The peudospin degree of 
freedom arises from the two-sublattice structure. Our analysis makes clear the mechanism how a phase transition 
occurs from a topological insulator to a band insulator under increasing electric field. We propose a method to 
determine the critical electric field with the aid of diamagnetism of silicene. Diamagnetism is tunable by the 
external electric field, and exhibits a singular behaviour at the critical electric field. Our result is important also 
from the viewpoint of cross correlation between electric field and magnetism. Our finding will be important for 
future electro-magnetic correlated devices. 



INTRODUCTION 

Silicene has recently been synthesizedlU-Q] and attracted 
much attention fsU?']. It is a monolayer of silicon atoms form- 
ing a two-dimensional honeycomb lattice. Almost every strik- 
ing property of graphene could be transferred to this innova- 
tive material. Indeed, its low-energy dynamics is described 
by the Dirac theory as in graphene. However, Dirac electrons 
are massive due to a relatively large spin-orbit (SO) gap of 
1.55meV in silicene, where the mass can be controlled by 
applying the electric field perpendicular to the silicene 
sheetldj]. A novel feature is that silicene is a topological 
insulatorl's'l, which is characterized by a full insulating gap 
in the bulk and helical gapless edges. 

Silicene undergoes a topological phase transition from a 
topological insulator to a band insulator as | increases and 
crosses the critical field E^, as has been shown|6] by exam- 
ining numerically the emergence of the helical zero energy 
modes in silicene nanoribbons. In this paper we present an 
analytic result by calculating the topological numbers based 
on the effective Dirac theory. We show that the origin of the 
topological numbers is a pseudospin meron in the momentum 
space. The pseudospin degree of freedom arises from the two- 
sublattice structure (Figll}- We also propose a simple method 
to determine experimentally the phase transition point with 
the use of the diamagnetism of silicene (FiglT]i- 

The magnetism of conventional metal is composed of the 
Pauli paramagnetism due to the spin magnetic moment and 
the Landau diamagnetism due to the orbital motion of elec- 
trons. The free electron system exhibits paramagnetism since 
the magnitude of the spin component is larger than the orbital 
component. Contrarily, the Landau diamagnetism overcomes 
the Pauli paramagnetism in a certain condensed matter sys- 
tem. An extreme case is provided by graphene [8^-|l2|, where 
the orbital susceptibility has a strong singularity due to the 
gapless character of Dirac electrons. 

We calculate the magnetic susceptibility of silicene as a 
function of the electric field Ez ■ We expect to have a strong 
singularity at the critical electric field, \Ez\ ~ E^-, since 
Dirac electrons become gapless at this point|6]. Indeed, we 
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FIG. 1 : (Color online) Illustration of an experimental setting to de- 
termine the critical electric field iSci with the use of the magnetic 
susceptibility of silicene. Silicene consists of a honeycomb lattice 
of silicon atoms with two sublattices made of A sites (red) and B 
sites (blue). The two sublattice planes are separated by a distance. A 
strong diamagnetism emerges between a magnet and a silicene sheet 
at the topological phase transition point — E„. 



show that it diverges at the critical electric field at zero tem- 
perature (T = 0), though the divergence is round off at fi- 
nite temperature. However, it is clearly observable as long 
as k^T < j^Aso- Our result is important also from the 
viewpoint of cross correlation between electric field and mag- 
netism. In general electric-field-controled magnetism is rather 
difficult compared to magnetic-field-controled electricity. On 
the other hand, the former is desirable since we can precisely 
control electric field. Our finding will be important for future 
electro-magnetic correlated devices. 



TIGHT BINDING MODEL 

Silicene consists of a honeycomb lattice of silicon atoms 
with two sublattices made of A sites and B sites. The two 
sublattices are separated by a distance, which we denote by 2£ 
with £ ~ 0.23A. The silicene system is described by the four- 
band second-nearest-neighbor tight binding modelira [Tsi [l4l] . 
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FIG. 2: The band gap As as a function of the electric field Ez- The 
gap is open for Ez ^ ±-Ecr, where silicene is an insulator. It has 
been shown|@] that it is a topological insulator for \Ez\ < E^ and 
a band insulator \Ez\ > E^- Thus there occurs a topological phase 
transition at |i5z| = Ea-. 
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FIG. 3: (Color online) One-dimensional energy bands for a silicene 
nanoribbon. (a) The bands crossing the gap represent edge states, 
demonstrating that it is a topological insulator, (b) All states are 
gapped, demonstrating that it is a band insulator. 



where cj^ creates an electron with spin polarization a at 
site i, and / {{i-.j)) run over all the nearest/next-nearest 
neighbor hopping sites. The first term represents the usual 
nearest-neighbor hopping with the transfer energy t = 1.6eV. 
The second term represents the effective SO coupling with 
Aso — 3.9meV, where cr ~ {a^, <Ty,az) is the Pauli matrix of 
spin, with uij = +1 if the next-nearest-neighboring hopping 
is anticlockwise and = —1 if it is clockwise with respect 
to the positive z axis. The third term represents the second 
Rashba SO coupling with Ar2 — 0.7meV associated with the 
next-nearest neighbor hopping term, where /i; = ±1 for the A 
(B) site, and dij — dij/ \dij\ with the vector dy connecting 
two sites i and j in the same sublattice. 

We take a silicene sheet on the xy-plane, and apply the elec- 
tric field Ez perpendicular to the plane. There appear two ad- 
ditional terms in the Hamiltonian, 



(2) 



Dirac theory around the A',, point as 



with 



hii = AsoO"z + aAR2 {kyGx — k^Uy) 



(3) 
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where Tq is the Pauli matrix of the sublattice pseudospin, 
v-p — -^cd is the Fermi velocity, and a ~ 3. 86 A is the lattice 
constant. 

The band gap is located at the K and K' points. At these 
points the energy is exactly given by 



(5) 



where s = ±1 is the spin-chirality. It is given by s = Szf] 
when the spin Sz is a good quantum number. The gap is given 
by 2 1 A, {E,) I with 



A, (Ez) = -sXso + \lE, + ].^{tE,f + \l, 



(6) 



where the first term represents the first Rashba SO coupling 
associated with the nearest neighbor hopping, which is in- 
duced by external electric field[15, 16]. It is proportional to 
the external electric field, Xr\{Ez) cc Ez, and becomes of 
the order of 10/ieV at E^, = Xso/i = 17meVA~^ The sec- 
ond term is the staggered sublattice potential term oc 2£Ez 
between silicon atoms at A sites and B sites. The total Hamil- 
tonian is given hy H = Hq + He- We note that the first 
Rashba SO coupling term (oc Ari) is missed in the previous 
analysis[6]. 



As \Ez \ increases, the gap decreases linearly since Ari oc Ez 
and vanishes at the critical point l-Bzl = E^- with 
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±17meV/A, 
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and then increases linearly (FigjS) The correction due to the 
first Rashba coupling is extremely small, (Ari/2Aso)^ — 



DIRAC THEORY 

Electronic states near the Fermi energy are tt orbitals resid- 
ing near the K and K' points at opposite comers of the hexag- 
onal Brillouin zone. We also call them the points with 
?7 = ±. The low-energy effective Hamiltonian is derived from 
the tight binding model H — Hq + He- It is described by the 



SPIN CHERN NUMBER 

We have shown in a previous paper[6] that silicene is a 
topological insulator for l-Bz] < E^ and it is a band insulator 
for \Ez\ < Eci- by examining numerically the emergence of 
the helical zero energy modes in a silicene nanoribbon (Figj3]|. 
In this section we present an analytic discussion by calculating 
the topological numbers based on the effective Dirac theory. 
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FIG. 4: Illustration of a meron in momentum space. The pseudospin 
texture is a meron for each spin in each valley, which yields 1/2 to 
the Pontryagin number. 



The topological quantum numbers are the Chem number C 
and the Z2 index. If the spin Sz is a good quantum number, 
the Z2 index is identical to the spin-Chern number Cs modulo 
2. They are defined when the state is gapped, and given by 



C = C+ + C_, = i(C+-C_), 



(8) 



where C± is the summation of the Berry curvature in the 
momentum space over all occupied states of electrons with 
Sz — ±1. They are well defined even if the spin is not a good 
quantum number[17, 18]. In the present model the spin is not 
a good quantum number because of spin mixing due to the 
Rashba couplings Ari and Ar2, and the resulting angular mo- 
mentum eigenstates are indexed by the spin chirality s — ±.1. 
A convenient way of calculating the Chern number C and the 
Z2 index is to use the formula (O to the system without the 
Rashba couplings and then adiabaticallyswitching on these 
couplings to recover the present svstem lU?! Iisll . 

When we set Ari = and Ar2 = 0, the Hamiltonian (O 
becomes block diagonal. For each spin Sz = ±1 and valley 
rj = ±1, it describes a two-band system in the form. 



H 



(9) 



where 



dx = rjhvpkx, dy = hvpky, dz = nio, (10) 
with the Dirac mass 

mo^ sXso+eEz. (11) 

The summation of the Berry curvature is reduced to the Pon- 
tryagin index in the two-band system!. 19l, 

where d = d / \ d\ is the unit vector which specifies the direc- 
tion of d. It is equal to the number of times the unit sphere is 
covered upon integrating over the Brillouin zone. 

It is convenient to use the cylindrical coordinate in the mo- 
mentum space, where 



dx ± idy — ^/l — a^(k)e^ 



rid 



<j{k) (13) 
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The pseudospin texture ( fT3T l describes a vortex-like meron in 
the momentum space, as shown in FigH] The Pontryagin in- 
dex ( fT2b yields 



^ J d^k eijdiadjO ^ %gn{mD) J da. 



4tt 

Hence we find 



(15) 



Cs = J2l(C+-C-) = ^i^^o-e\Ez\), (16) 



where 8 is the step function, i.e., 

1 for \eEz\ < Aso 



for \iE,\> Aso 



(17) 



We have verified that the system is a topological insulator for 
\iEz\ < Aso and a band insulator for \iEz\ > Aso in the 
system without the Rashba interactions. The property remains 
true when they are switched on adiabatically. 



DIAMAGNETISM 

We proceed to discuss a possible experimental method 
to detect the phase transition point by measuring the mag- 
netic susceptibility. We apply homogeneous magnetic field 
B = V X A = (0, 0, -B) with B > along the z axis to 
silicene|7]. By making the minimal substitution, the Hamilto- 
nian is given by 



Hr, =Vp {riPxTx + PyTy) + T]Tzhll + lEzTz 



(18) 



with the covariant momentum Pi = fiki + eAi. We introduce 
a pair of Landau-level ladder operators. 



iB{Px+tPy) 

V2h 



V2h 



(19) 



satisfying [d,d^ = 1, where £b = ^/h/eB is the magnetic 
length. In the basis {-tpAt^ V'sf ^Ai, "05^}*' the Hamiltonian 



reads 



+ 
huj^d 



fujJcd^ 
-A^ {Ez) 
zAri 



-zAri 

A*! [Ez) 
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(20) 
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FIG. 5: Susceptibility x as a function of electric field Ez for vari- 
ous temperature keT/Xso ~ 0.01, 0.1, 1, 10. It has a sharp peak at 
\Ez\ = £cr for /cbT/Aso < 0.1, and diverges asT -^0. 



at the K point, with Wc = V^hvp/^B- Here the diagonal ele- 
ments A J. (Ez) are 

Al{E,) = ±Xso + £E,t, (21) 

with the sublattice pseudospin t^- We note that A± (Ez) — 
A% {Ez)ifv/e setARi =0. 

The magnetic susceptibility is defined by 

Y = lim (22) 

where Af is the magnetization and B is the external magnetic 
field. The general formula for the orbital magnetic suscepti- 
bility of Bloch electrons is given by||20j], 



-gv 



^ ^ Tr {Gv,GvyGv,Gvy) , (23) 



n k 



where Vi = dH/dki, /3 = l/fceT, and G is the temperature 
Green function 



G {k,LUn) = {ihuJn - H) \ 



(24) 



with Hun = {2n + 1) tt//? being the Matsubara frequency. 

By making the Taylor expansion of (|23] | with respect to Ari 
and Ar2, the Rashba terms are found to yield second order 
corrections, Sx — o(A^j/?i^i;p) + o(A^2/^^''^f)' hence we 
neglect the Rashba terms in what follows. In this approxima- 
tion the spin Sz is a good quantum number 

Integrating out the wave number k of the matrix trace of 
), we have 



Y,^r (Gv^GvyGv^^Gvy) = ^ 
k s=± 

Using the formula of the infinite sum, 

1 1 



(25) 



E 



A? 2A,fcBT 



■ tanh 



2fcBT' 



(26) 



we obtain the magnetic susceptibility at finite temperature 

e'^ffv^ ^ 1 .A 



(27) 



where gy = 2 represents the valley degree of freedom. We 
show X {T, Ez) as a function of Ez for typical values of T. 



We can make a precise determination of the critical elec- 
tric field Ea- based on this formula. Let us consider, for ex- 
ample, a situation where a magnet is placed parallel to a sil- 
icene sheet (FigHli. The magnetization of silicene is given by 
M = — where B is the external magnetic field made 
by the magnet. The magnet feels a strong repulsive force 
at Ez = Ecr, since we have As — > and M — > — oo as 
\Ez\ -> E'er at r = 0. Even for finite temperature, x {T, Ez) 
has a sharp peak at \Ez \ = E^- for k^T < Aso/10 as in Fig|5] 
Such a strong repulsive force can be detected mechanically. 
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